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In the one-way model of quantum computing, quantum algorithms are implemented using 
only measurements on an entangled initial state. Much of the hard work is done up-front 
when creating this universal resource, known as a cluster state, on which the measurements 
are made. Here we detail a new proposal for a scalable method of creating cluster states 
using only a single multimode optical parametric oscillator (OPO). The method generates 
a continuous-variable cluster state that is universal for quantum computation and encoded 
in the quadratures of the optical frequency comb of the OPO. This work expands on the 
presentation in Phys. Rev. Lett. 101, 130501 (2008). 



I. INTRODUCTION 

The invention of one-way quantum computing (QC) introduced a new paradigm for quan- 
tum information processing [1], whereby the unitary evolution required in the traditional circuit 
model [2] is replaced by measurements on a highly entangled cluster state [3]. In this model, 
quantum information exists virtually in the cluster state and is manipulated through adaptive 
single-qubit projective measurements. The measurement backreaction on the entangled state ef- 
fectively induces dynamics, and hence quantum computation. In this way, the cluster state acts 
as a quantum field-programmable gate array, or a quantum "substrate" [1], on which a quantum 
circuit is inscribed and simultaneously evaluated via measurements. 

Because it eliminates the issue of dynamic control, one-way QC is an especially attractive model 
for certain physical systems where creating cluster states is natural. Natural, in this context, 
means that large cluster states (say, with more than 100 qubits) can be created in just a single 
step, without the need to build the cluster piecewise before the computation begins. The most 
prominent example of such a system is an optical lattice with cold atoms [4, 5], but other proposals 
exist, e.g., using solid-state qubits [6]. These systems provide ambitious "top-down" approaches to 
quantum computing, where the issue of scaling the number of qubits is addressed up front. This 
is in contrast to the more widely attempted "bottom-up" approaches which initially emphasize 
one- and two-qubit control [7, 8, 9], such as ion traps [10], linear optics [11], and solid state 
implementations [12, 13], among others [14]. 

In this paper, we provide a detailed proposal for another candidate for creating large cluster 
states efficiently: a single optical parametric oscillator (OPO) driven by a multifrequency pump. 
This result was first reported in Ref. [15]. This article is an exhaustive description that fills in 
details of the derivations and proposed implementations. 

The eigenmodes of an optical cavity (e.g. two mirrors facing each other), whose frequencies 
satisfy a constructive interference requirement, define an optical frequency comb (OFC), so called 
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because the frequencies are evenly spaced. The OFC has outstanding classical coherence properties, 
as illustrated by the revolutionary and now ubiquitous use of million-frequency, mode- and carrier- 
envelope-phase-locked femtosecond lasers in time and frequency metrology [16, 17]. When the 
electromagnetic field is quantized, this comb of eigenmodes forms a set of independent harmonic 
oscillators which we propose to use as the core of a quantum computer by entangling them into a 
cluster state. The goal being to generate cluster-state entanglement between these quantum modes, 
or qumodes [67], we first augment the optical cavity with a nondissipative, second-order nonlinear 
gain medium, thus converting it into an optical parametric oscillator (OPO). With an appropriate 
nonlinear medium and a pump beam with the right frequency content, the comb of qumodes 
will encode a large square lattice continuous- variable cluster state [15]. Remarkably, this can be 
achieved using just a single OPO — and with a constant number of pump frequencies, regardless 
of the size of the cluster. Compared to existing proposals [18, 19], these modest experimental 
requirements, together with the naturally large set of qumodes inside the comb mean that this 
approach has strong potential for scalability. Detailing this construction is our main goal. 

The paper is organized as follows. In Section II, we discuss the basic physics of the OPO. In 
Section III we define continuous- variable (CV) cluster states, drawing distinctions with their qubit 
counterparts, and discuss previous literature dealing with the creation of CV cluster states. In 
Section IV we outline the program of creating a CV cluster state from just a single OPO and a 
multifrequency pump [20] and introduce a simplifying ansatz [21]. We demonstrate the power of 
this approach in Section V, where we detail the construction of a ring-shaped CV cluster state that 
is universal for single- mode transformations. The ring construction contains all of the essential 
ideas for the main result in Section VI, which was first discussed in Ref. [15]: a fully universal 
CV cluster state from a single OPO using only 15 pump frequencies, regardless of the size of the 
cluster. We detail a proposed experimental implementation of this scheme in Section VII. The 
important but unresolved issues of error correction and fault tolerance are discussed in Section IX, 
and we conclude in Section X. 



When an OPO is driven by a pump field with frequency and wave vector (u> p , k p ), the nonlinear 
medium inside promotes the up- or downconversion of pump photons into pairs of signal photons 
[(u m , k m );(w n , k n )], if the phasematching conditions 



hold, in particular so long as the photons have frequencies within the phasematching bandwidth 
of the medium. This interaction is also dependent on other optical quantum numbers, such as the 
polarization and the transverse spatial mode. We assume for now that only one polarization per 
signal frequency gets phasematched for a given pump frequency; we will relax this restriction later. 
The phasematching is assumed collinear. 

It is by now well established that such interactions lead to bipartite CV entanglement between 
the signal qumodes [22] both below [23] and above [24, 25, 26] the OPO threshold. Any single such 
conversion, a.k.a. parametric, process has an interaction-picture Hamiltonian of the form 7i = 
ihn(alna,n — a m a n ), where k = fix is a nonlinear coupling strength (squeezing parameter per 
unit time), with x the nonlinear susceptibility of the medium and /3 the complex amplitude of 
the pump (assumed undepleted and classical), and a„ is the creation operator for photons with 
frequency uj n . We define the amplitude and phase quadrature operators for the k th qumode as 
Ik = (®k + ^fc)A/2 and pk = i{a\ — d fe )/-\/2. The sign of k is set by the pump phase equal to 0, 



II. OPTICAL PARAMETRIC OSCILLATORS 




(1) 
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in which case the Hamiltonian Ti generates quantum squeezing of the phase-sum and amplitude- 
difference simultaneously, which constitutes a realization of the Einstein-Podolsky-Rosen paradox 
[27, 28] and therefore bipartite entanglement. Similarly, if the pump phase is tt, there will be 
simultaneous amplitude-sum and phase-difference squeezing. Note also that these two possibilities 
coincide respectively with parametric amplification and deamplification of the signal fields when 
they are seeded by coherent states. 

A single pump frequency can simultaneously generate multiple pairwise interactions, and there- 
fore multipartite entanglement, if there are multiple pairs of signal frequencies that satisfy Eq. (1) 
within the phasematching bandwidth of the crystal [21, 29, 30]. In this case, we can write the full 
interaction-picture Hamiltonian as 

H(G) = ihK Gmn{a\ n ai - a m a n ) . (2) 

m,n 

This Hamiltonian is written explicitly as a function of a matrix G (with entries G mn ), which is the 
adjacency matrix describing the pairwise coupling strength between each pair of signal modes in 
units of k. The matrix G is real and symmetric, with positive (negative) entries corresponding to 
down- (up-)conversion. 

So far, this is a fairly general Hamiltonian, since we have said nothing about the ordering of the 
signal modes. In an OPO the signal modes are equally spaced in frequency, forming what's called 
an optical frequency comb, so we can choose to order them sequentially like this: 

oj n = O + nAw , (3) 

where £1 is an offset frequency (discussed below), and Au; is the free spectral range of the OPO 
(i.e., the mode spacing). We neglect the effects of optical dispersion in this paper. Not all values 
of n or m correspond to entries within G, however, since the only modes involved in the interaction 
are those within the phasematching bandwidth. We therefore choose Q such that UJ\ corresponds 
to the lowest-frequency mode within the phasematching bandwidth and lon to the highest, which 
allows us to label the rows and columns of G each from 1 to N (inclusive), corresponding exactly 
to the cavity modes available for interaction. With this convention, the phasematching condition 
of Eq. (1) becomes n + m = (const.), which selects a skew-diagonal band to be the only nonzero 
entries in G. (The particular band selected is determined by the relationship of the pump frequency 
to the phasematching bandwidth.) The entries in this band are all equal and given in units of k, 
with the sign indicating downconversion if positive or upconversion if negative, as above. 

This construction generalizes easily to the case where the OPO is pumped with multiple frequen- 
cies simultaneously. In this case, each pump frequency corresponds to a different phasematching 
condition, Eq. (1). Thus, each corresponds to a different skew-diagonal band in G. Since we want 
all such interactions summed together as in Eq. (2), this leads us to consider exclusively those G's 
that are of Hankel form — i.e. those that have constant skew-diagonals. For convenience later, we 
introduce a shorthand notation for Hankel matrices of any size. A Hankel matrix is completely 
specified by its first row and last column, so we will collect these terms into a vector, setting aside 
by slashes the middle element common to both. The shorthand form of a given Hankel matrix is 
unique and represents the entire matrix itself. As such, the following is a literal equality: 

(a b c\ 
ted . (4) 
c d ej 

While this shorthand form can be used to represent any Hankel matrix, it has an additional 
advantage when used to represent G in that it corresponds simply and directly to the spectrum of 
the generating pump beam, as illustrated in Figure 1. 
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FIG. 1: Hankcl shorthand and pump specification. A Hankcl matrix is uniquely specified by the entries 
along the top and down the right side. We collect these entries into a shorthand vector, representing the 
entire matrix itself, with the top-right entry set off with slashes. When G is Hankel, its shorthand vector 
immediately specifics the pump spectrum required to implement Ti(G). Each nonzero entry in the shorthand 
vector denotes the amplitude of a frequency in the pump, each of which generates CV entanglement between 
pairs of qumodes in the OFC symmetric about half that frequency. This accounts for all couplings prescribed 
by G. 



Polarization and Block-Hankel Matrices 

The description above assumes that only one interaction per pump frequency is phasematched 
in the crystal. Appropriate additional quantum numbers (polarization, spatial mode, etc.) for 
the pump and signal modes are implied. An important extension of this scheme is when there 
are multiple modes that are frequency degenerate and yet phasematched independently by the 
crystal. The simplest example of this is when the crystal phasematches multiple polarization- 
dependent interactions, but other possibilities include wavevectors and transverse spatial modes. 
The simultaneous quasiphasematching [31, 32] of four polarization-sensitive interactions ZZZ, ZYY, 
and YZY/YYZ (where the first letter is pump polarization) has been demonstrated experimentally 
in a periodically poled KTiOP0 4 (PPKTP) crystal [33]. 

In such a case, the Hankel restriction on G can be relaxed to block-Hankel, with the block 
size given by the number of frequency degenerate modes. Take the example of polarization- 
dependent interactions discussed above. If M frequencies are phasematched by the crystal, with two 
polarization-modes each, then there are N = 2M total modes, ordered ly, lz> 2y, 2z, • • • , My, Mz, 
with Y' and 'Z' labeling the two polarizations. Extensions including more quantum numbers are 
possible for setups that involve more complicated frequency-degenerate phasematching. The im- 
portant thing is that each frequency label is completed by a set of additional mode labels, with 
the same ordering of additional quantum numbers within each set. 

Returning to the simplest case of polarization-dependent phasematching (which is all we will use 
in this paper), a single pump frequency can now generate several interactions at once, depending 
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FIG. 2: Polarization-sensitive interactions and block-Hankel matrices. This is a generalization of Fig- 
ure 1 to include polarization-sensitive interactions. Each of the four frequencies now has two independent 
polarization- modes that can be coupled within the crystal simultaneously. The Hamiltonian adjacency ma- 
trix G is now 8x8 and is no longer strictly of Hankel form but rather is block-Hankel with 2x2 blocks. These 
blocks along the top and right side can still be collected into a shorthand vector, as shown, where p = ( £ £ ) , 
a = ( § c ) ' an< ^ represents a block of zeros. The nonzero blocks each specify the amplitude and polarization 
angle in the (YZ)-plane of a single pump frequency when used with a PPKTP crystal [33]. 



on its polarization. Using the crystal discussed above [33], a monochromatic pump implements a 
skew-diagonal band of 2 x 2 blocks in G. The particular form of the block is determined by the 
pump's polarization angle in the (YZ)-plane. Multiple polarized frequencies (possibly at different 
angles) in a single pump will produce a block-Hankel G, again with 2x2 blocks, as illustrated in 
Figure 2. Furthermore, the Hankel shorthand discussed above may be extended intuitively to block- 
Hankel matrices. The shorthand vector is now comprised of blocks instead of single entries and, 
when applied to G, still directly specifies the pump spectrum, including the necessary polarization 
of each frequency in the beam. 



III. CONTINUOUS- VARIABLE CLUSTER STATES 



The relationship between the qubit-based and CV-based formalisms is by now well-established, 
and spans the areas of QC [34, 35], quantum error correction [36, 37], cluster states [38], and one- 
way QC [18]. Given a collection of N qumodes, we collect the single-mode quadrature operators 
into column vectors given by q = (q±, . . . , qN) T and p = (pi, . . . ,Pn) T ■ Then the definition of a 
CV cluster state [18, 38, 39] is any Gaussian state whose quadratures satisfy, in the limit of large 
squeezing, 



p(r) — Aq(r) r ~* 00 > , 



(5) 



where r > is a squeezing parameter and A is the adjacency matrix of the CV cluster's graph. 
Because they are equivalent descriptions, we can speak interchangeably about the graph of the 
cluster state and its adjacency matrix A. The usual choice of A for universal cluster-state QC is a 
square-lattice graph, but others are possible as well. Many regular lattices [40] and also irregular 
lattices above a percolation threshold [41, 42], and even certain more exotic states [43] are now 
known to be universal resource states. 

The infinite-squeezing limit of Eq. (5) is an important difference from the standard (discrete- 
variable) cluster state; it is not achievable by any finite-energy state, so we must be content with 
finitely (albeit strongly) squeezed states as an approximation. This means that errors due to finite 
squeezing are inherent in CV cluster state QC, a point we return to in Section IX. 

Previous proposals for creating CV cluster states called for O(N) vacuum squeezers and O(N) 
quantum non-demolition gates, which can be implemented by experimentally arduous inline squeez- 
ers (e.g., seeded OPOs) [44, 45], although realization with more convenient offline, vacuum squeezers 
was recently achieved [46]. An elegant simplification is provided by the Bloch- Messiah decompo- 
sition [19, 45], which reduces the experimental requirements to N OPOs followed by a 2V-port 
interferometer. This method has already been used successfully to demonstrate the creation and 
use of four-mode CV cluster states for simple quantum information processing tasks [47, 48, 49]. 
In contrast, our method requires a single OPO and no interferometer. As we shall see, the pump- 
ing field's complexity is a constant with N . The required nonlinear medium is sophisticated yet 
feasible with current technology and has been demonstrated experimentally [33]. 

Quantum computation on the CV cluster state [18] proceeds analogously to the case of qubit 
cluster states [1]. Adaptive measurements are performed on the cluster to induce computation. 
Homodyne detection alone is sufficient to implement all Gaussian operations, which are the analogs 
of Clifford gates for qubits. In addition, the required basis for a given homodyne measurement does 
not depend on the outcomes of other homodyne measurements. Thus, measurements implementing 
Gaussian operations can be performed in any order, and even in parallel. The availability of a non- 
Gaussian measurement (e.g., photon counting) is essential for universal QC [18]. Non-Gaussian 
measurements, like non-Clifford gates for qubits, are not parallelizable in general because future 
measurement bases will depend on the outcome of this measurement. 

IV. CV CLUSTER STATES FROM A SINGLE OPO 

In this section, we review the relationship between the coupling matrix G of the OPO Hamil- 
tonian in Eq. (2), and the adjacency matrix A of the created CV cluster state in Eq. (5). Our 
treatment differs from the one in Ref. [20] in that we make an additional simplifying ansatz that 
was first made explicit in Ref. [21]. This analysis also makes explicit the scaling with respect to the 
squeezing parameter in Eq. (5). A more general analysis of the underlying mathematical structures 
can be found in Ref. [50]. 

As we have seen in Section II based on physical arguments, we should specialize to G's that 
are Hankel or block-Hankel, as in Figures 1 and 2. We make the simplifying ansatz that there are 
no degenerate squeezing interactions within the OPO, which, with the (block-)Hankel requirement 
ensures that G is the adjacency matrix of a bipartite graph. Any bipartite G can be factored into 
a tensor product as 

G = A ® X (6) 
where X = ( ? q) , an< ^ is a block-Hankel matrix having the same size blocks as those of G. 
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For example, consider the configuration depicted in Fig. 1. We can factor G as 



[0, 0, 0/a/0, b, 0] 



/0 a\ 
a 
a 6 

\o 6 0/ 



a 
a b 



1 

1 



[0/a/b] ® [0/1/0] = A (8) X , (7) 



where we move freely between the matrix form of G and the Hankel shorthand discussed above. 
A Hankel G with an even number of modes has such a factorization if and only if it contains no 
single-mode squeezing — or equivalently, if it has only zeros on the main diagonal. 

We now introduce the additional ansatz of Ref. [21] and assume that G T G = G 2 = 1, or 
equivalently, that Ag = 1. (Here and throughout, the size of the identity matrix represented by the 
symbol 1 is given by the context. The same holds for the zero matrix denoted by 0.) By assuming 
that G is orthogonal, we greatly simplify its relationship to A. In fact, we will simply define 



A := X <g> A r 



(8) 



and show that it is the adjacency matrix for the CV cluster state generated by TC(G) [20]. This 
definition would imply that A and G are isomorphic graphs, since reversing the tensor product order 
of an adjacency matrix corresponds merely to a particular permutation of node labels. Continuing 
with the example from above, 



X<g> A 



o 



A 
A 



/0 a\ 
a b 
a 

\a b 0/ 



fO a\ 
a 
a b 

\a b 0/ 



X 



(9) 



We denote equivalence up to permutation by =, which in this case involves exchanging the labels 
of modes 2 and 3. 

To show this correspondence, we begin by collecting the quadrature operators q and p into a 
single concatenated vector v = (q±, . . . ,qN,Pi, ■ ■ ■ ,Pn) T which represents the vacuum (i.e., their 
measurement statistics are those of the vacuum state). In the Heisenberg picture, symplectic linear 
transformations on this vector correspond to Gaussian transformations in Hilbert space [51]. The 
unitary operation exp(-itTt) of the OPO Hamiltonian from Eq. (2) acts on the qudrature operators 
(in the Heisenberg picture) through the following linear symplectic matrix acting on v: 

orC 



U 









-rC 



(10) 



where r = Kb. We additionally allow permutation of the qumode labels and phase-shifting of in- 
dividual qumodes; these are completely passive operations and don't change the required physical 
setup used to generate the state. Using the permutation freedom, we choose to reorder the quadra- 
tures so as to reverse the tensor product order, as in Eq. (8), effectively converting G to A, and 
resulting in the redefinition 



U 







(11) 



Using G 2 = A 2 = 1, we can evaluate U using the identity e ±rA = cosh(r)l ± sinh(r)A, resulting in 



U 



/cl 
sA 


V o 



sA 
cl 







cl 

-sA 






-sA 
cl 



(12) 
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where c = cosh(r), s = sinh(r), and each block is the same size as Ao- We now phase shift half of 
the qumodes by —ir/2 using the symplectic matrix 



T 



/l \ 

-1 

10 

\0 1 / 



(13) 



This operation does not need to be implemented actively since it is simply a redefinition of the 
quadratures (q — ► — p and p — > q) for the phase-shifted modes. This results in 



TU 



/ cl sA 


\sA cl 




sA 
cl 




\ 

-cl 

-sA 
/ 



(14) 



To verify that Eq. (5) holds, it is sufficient to check that when we act on the transformation 
matrix TU with the rectangular matrix 



(-A 1) 





-A 



-A 1 
1 



we get a matrix that vanishes as r — > oo. We find that, in fact, 

1 A \ r->c 



-A 1)TU = e~ r 



A 1 







(15) 



(16) 



This shows that the CV cluster-state relation, Eq. (5), is satisfied. 

We have shown that if a Hamiltonian coupling matrix G can be found that satisfies the (rather 
restrictive) assumptions above — namely, that it be (block-)Hankel, orthogonal, and bipartite — then 
it creates a CV cluster state with an identical adjacency matrix, up to permutation of labels and 
local redefinition of quadratures. In the next two sections, we will provide examples of such G's 
that yield CV cluster states with very interesting graphs using only a constant number of pump 
frequencies. In Section V, we show how to construct a ring graph which is universal for single mode 
transformations, and in Section VI we show how to construct a torus, which is fully universal for 
quantum computation. In what follows, since they are equivalent, the distinction between G and A 
will be eliminated and A used to represent both. 



V. SINGLE-MODE-UNIVERSAL CV CLUSTER STATE 



We first illustrate the main ideas on a simpler graph that is universal for single-mode operations 
(but not for universal QC, which requires multimode operations). We desire such a CV cluster state 
that also has a bicolorable graph [20] and whose adjacency matrix is orthogonal and Hankel for 
experimental simplicity [21]. Orthogonality of the adjacency matrix (AA T = 1) for an undirected 
graph (A = k T ) yields A 2 = 1, or 

(A 2 )jk = ^A j iA ik = 5 jk . (17) 
l 

Eq. (17) has a geometric interpretation: (A n )jk represents the sum of the weights of all n-length 
paths from node j to node k, where the weight of such an "ra-path" equals the product of all edge 
weights along the path. Eq. (17) enforces two conditions: 
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FIG. 3: Interpretation of a matrix-weighted edge. A matrix-weighted edge connecting two macronodes 
specifies the graph connecting the underlying physical nodes. In this example, the macronodes are connected 
by a 2 x 2 matrix-valued weight (left). Each macronode itself therefore contains two physical nodes, and 
the real- valued weights connecting them are specified by entries in the matrix (right). 



(1) all 2-paths that begin and end on the same node have weights that sum to 1, and 

(2) all 2-paths that link different nodes have weights that cancel out. 

The simplest graph to try would be a line graph [1] but such a graph necessarily fails condition (1) 
because of irregularity at the boundaries. The next natural graph to consider is therefore a ring 
graph, which is a regular graph, but orthogonality is still prohibited because each node is connected 
to its next-nearest neighbor by only one 2-path [68], for which the sum in Eq. (17) collapses to a 
single term (which must be nonzero), thus violating condition (2). 

While real-valued weights cannot satisfy Eq. (17) for a ring graph, matrix-valued weights can, 
as shown below. In such a case, Eq. (17) becomes 

J2Ni*ik = 8jkl, (18) 
i 

where the "entries" Aji are themselves mx m matrices. This means A is now an adjacency matrix 
on an (miV)-node graph. On the other hand, treating the m x m blocks as single entries, A is the 
matrix- weighted adjacency matrix for what we call a supergraph, which has macronodes consisting 
of m individual nodes each. Figure 3 illustrates the meaning of matrix-weighted edges between 
two macronodes. 

Promoting the ring graph to a supergraph with 2 nodes per macronode (i.e., m = 2), we can 
choose our weighting to consist of two orthogonal projectors over M 2 , labeled w^, with A,-j-_|_i = tt + 
for j even and Ajj+i = tt~ for j odd, where the total number of macronodes N is even, and 
index addition is modulo N. By reasons of symmetry and connectivity of the underlying graph, 
we choose tt^ to be 

- ---*(-;)• (i9) 

with ± standing for ±1. With these weights, the supergraph is still a ring, but the actual graph 
on physical nodes has a more complicated "crown" structure, as shown in Figure 4. Either by 
direct verification of Eq. (18) or by noting that the geometric conditions (1) and (2) are satisfied, 
this weighting of the ring supergraph results in an orthogonal A. Given a CV cluster state with 
this graph, one can perform local measurements on one of the two rings of physical nodes and 
measure it down to a CV cluster state [18] having a simple ring graph with equal weights, which is 
universal for single-mode operations using one-way QC methods [1, 39]. The only thing remaining 
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FIG. 4: Matrix- valued weights and supergraphs. The matrix- valued weights ■k ± , defined in Eq. (19), connect 
the macronodes of a "ring" supergraph (left). The entries in ^ specify the real- valued weights in the actual 
"crown" graph (right) that connects the underlying physical nodes. Measuring q for each of the physical 
nodes in the top layer of the crown leaves the bottom layer in a uniformly weighted ring-graph CV cluster 
state. 

is to show that the adjacency matrix for the crown can be put into Hankel or block-Hankel form 
by an appropriate numbering of the nodes. 

It is easiest to begin at the supergraph level. The ring is a circulant supergraph, which means 
it can be represented by a circulant adjacency matrix simply by numbering the nodes sequentially 
around the ring. What we want, however, is the skew-circulant form, a special case of Hankel 
where the skew-diagonal bands "wrap around" at the edges of the matrix. We can renumber the 
nodes to achieve this form as follows. Starting with an arbitrary node, label it 1. Move two nodes 
clockwise, and label that node 2. Continue labeling every other node sequentially in a clockwise 
fashion until you reach node N/2. Then, pick one of the remaining nodes to label as N/2 + 1, 
move two nodes counterclockwise, and label that node N/2 + 2, continuing in this fashion until 
all nodes are labeled. The result is a skew-circulant adjacency matrix at the supergraph level, 
with matrix- valued weights given by Eq. (19). Such a matrix is block-skew-circulant (and thus 
block-Hankel) at the physical node level: 

A = [0, . . . , 0, 7r + , 0/ir~ / 0, . . . , 0, 7r + , 0] , (20) 

N-3 N-3 

where the zeroes are 2x2 blocks of all zeros, and = indicates equality up to renumbering of nodes. 

We could stop here, since we have an adjacency matrix that is block-Hankel with 2x2 blocks. 
The Hamiltonian TC(A) from Eq. (2) can be implemented using a pump beam with three polarized 
frequencies and a PPKTP crystal [33] using the method described in Section II. We will find it 
useful, however, to introduce a renumbering scheme that can be used to make the crown fully 
Hankel because a similar scheme will be used in the next section to reduce a block-Hankel matrix 
with 4x4 blocks to one with 2x2 blocks. 

Since A is an N x N matrix of 2 x 2 blocks, any index j can be thought of as an ordered pair 
[m, c] := 2m + c = j, with m ranging from to N — 1 and c being either or 1. The permutation 
that reverses the order of the subindices, [m, c] — > [c, m] := Nc + m = j' , can be used to convert A 
into fully Hankel form: 

A = ^[O^ -3 , 1, 0, 1, O^ -3 , 1, 0/—1/0 N ~ 3 , 1, 0, 1, O^ -3 , 1, 0] , (21) 
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where O n represents a string of n zeros. Being fully Hankel, Eq. (21) can be used to directly read 
off the seven pump frequencies needed to implement this CV cluster state in a single 0P0 using an 
optical frequency comb of qumodes by the correspondence shown in Figure 1. Notice that since the 
only effect of increasing N is to add more zeros to this Hankel vector — shifting, but not increasing 
in number, the required pump frequencies — the complexity of the pump is constant with respect 
to N, making this construction extremely scalable. 

The CV cluster state described above therefore fulfills the four criteria desired in this section: 
it is bipartite, it has an orthogonal adjacency-matrix representation, that representation is also 
Hankel, and it can be used to implement any single-mode operation. The first property satisfies 
the necessary condition for generation in a single OPO [20]. The second satisfies our desire for 
mathematical simplicity, allowing us to use the cluster-state graph A without modification as the Ti- 
graph from Eq. (2). The third guarantees simple experimental implementation [21], with constant 
pump beam complexity. The fourth says that single-mode measurements can be used to effect any 
single-mode unitary operation. While not universal for multimode QC over continuous variables, 
the construction described here has introduced all the essential concepts for a construction that 
achieves this goal, described next. 



The construction in the previous section gave a spectacularly efficient method of generating a 
CV cluster graph with one-dimensional topology (such as a line or ring), hence only suitable for 
single-mode operations. In universal one-way QC, however, operations involving more than one 
qubit/qumode require graph connections between such one-dimensional "quantum wires." The 
natural candidate for a fully universal CV cluster-state graph is the square lattice — the original 
graph proposed for qubit cluster states [1] — but, similar to the case made against the line graph, 
a square lattice does not admit an orthogonal adjacency matrix because it is not regular at the 
boundaries. A logical response is then to impose toroidal boundary conditions, linking nodes on 
one side to those on the opposite side. However, the resulting toroidal square lattice suffers from 
the same problem as the simple ring: there exist pairs of nodes connected by exactly one 2-path, 
prohibiting an orthogonal representation. 

We therefore take the same measures and promote the toroidal lattice to a supergraph. Recall 
that orthogonality requires that A 2 = 1, equivalent to the aforementioned two geometrical condi- 
tions (1) and (2), which are local on the (super-)graph and thus are not affected by the overall 
topology. In this case, the supergraph locally looks like a square lattice, which has degree four, 
and Figure 5 illustrates that four mutually orthogonal projectors are sufficient to satisfy the geo- 
metric conditions. Once again motivated by symmetry and connectivity of the underlying graph, 
we choose the weights to be 4 x 4 projectors constructed from the 2x2 projectors in Eq. (19): 



We now possess an orthogonal supergraph, with 4-node macronodes, that looks locally like a square 
lattice but has toroidal topology. The final task is to convert the resulting adjacency matrix into 
Hankel form. 

Some freedom remains in the way that the lattice is rolled up into a torus that will turn out 
to make a Hankel form of the supergraph particularly easy to implement. The freedom lies in the 
fact that we can twist the lattice as we roll it up: instead of connecting the last node of each line 
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FIG. 5: Four-color solution to the geometric orthogonality conditions. Labeling the edges of a large lattice 
with four colors as shown, where each color represents one of four orthogonal projectors over R 4 , satisfies the 
geometric orthogonality requirements (see text). From each node protrudes exactly one edge of each color, 
satisfying condition (1), and any 2-path connecting distinct nodes (several shown) traverses edges with two 
different colors, which has exactly zero weight, satisfying condition (2). 



of the lattice to the first node in that same line, we can instead connect it to the first node of a 
different line — one that is n lines away from the original one — generating a twist of n lines in the 
torus. Finding a Hankel form for the ring relied on the fact that the graph is circulant, along with 
a suitable node renumbering procedure, to generate a skew-circulant adjacency matrix. While we 
could start with a toroidal lattice and try to find a twist that admits a Hankel form for it, we will 
instead explicitly construct a circulant graph and show that it is in fact a twisted toroidal lattice. 
We will construct the graph from two different circulant "threadings" of M 2 nodes, corresponding 
to the two different dimensions of the lattice. 

The simplest circulant threading is just an M 2 -node ring graph (we require M to be even). 
Alternating between two of the projectors from Eq. (22), say IT 2 and IT 3 , this is exactly the ring 
supergraph from the previous section (except now with 4-node macronodes). We then apply our 
second "threading" by additionally connecting each node in the ring to the two nodes that are M+l 
steps away, alternating with the remaining two projectors, in this case IT and IT 1 (see Figure 6). 
Since M+l and M 2 are relatively prime, this procedure creates another complete cycle through the 
ring, this time traversing every node by taking steps of size M + l. The result is a circulant graph 
where each node now has four neighbors. While tracing any path through the nodes, each step 
will move in one of the two threading "directions," moving by ±1 or ±(M + 1). These correspond 
to moving horizontally or vertically, respectively, through the lattice. The circulant nature of the 
graph guarantees toroidal boundary conditions. Thus, we are left with a twisted toroidal square 
lattice supergraph, shown in Figure 7. 

Since this graph is circulant, we can use the ring renumbering trick from the previous section 
(with N — > M 2 ) to generate a skew-circulant adjacency matrix for this supergraph, which is block- 
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FIG. 6: Circulant embedding of a twisted toroidal lattice. The 16 white spheres each represent a four-node 
macronode, while the colored edges each represent a matrix- valued weight IP from Eq. (22), with II red, 
II 1 yellow, II 2 blue, and II 3 green. The toroidal "axes" are identified with moving along the "red-yellow 
direction" and the "blue-green direction." The circulant construction includes a one-unit twist in each 
toroidal dimension and generalizes easily to more nodes (see text). 



skew-circulant (and thus block-Hankel) with 4x4 blocks at the level of the physical nodes: 

[o ll ,n 1 ,o^,n°,o il ,n 3 ,o/n 2 /o u ,n 1 ,o^,n°,o ti ,n 3 ,o] , (23) 

where u = (M — 1) and v = (M 2 - 2M - 3). Once again, the resulting matrix is block-Hankel 
instead of fully Hankel, this time with 4x4 blocks. For the ring, we used a renumbering of the 
nodes to convert the 2 x 2-block-Hankel matrix from Eq. (20) into the fully Hankel matrix in 
Eq. (21). An equivalent renumbering can be done in this case to convert this 4 x 4-block matrix 
to one with 2x2 blocks, which almost preserves the Hankel structure (we'll show how to fix it 
below): simply treat the 4x4 blocks of A as 2 x 2-block matrices of 2 x 2 blocks, and apply the 
renumbering from the previous section (this time, with N — > 2M 2 ). The resulting matrix is not 
completely block-Hankel, but it is close. We can fix this by defining A' as in Eq. (23) but with 
the sign nipped in the next-to-last entry (n 3 — ► — n 3 ). While A is block-skew-circulant, A' is only 
block-Hankel (but still orthogonal): 

a' = [o u , n 1 , o u , n°, o", n 3 , o/n 2 /o u , n 1 , cr, n°, o", -n 3 , o] , (24) 

The important thing, though, is that if we apply the renumbering procedure to A', the result is a 
2 x 2-block-Hankel adjacency matrix: 

A' ^ [0 S , 7T~, 0*, 7T+, S , 7T+, 0, 7T~, S , 7T", 0*, 7T+, 0* , 7T~, 0/vr+/ 

s ,7r",0 < ,7r + ,0 s ,7r + ,0,7r",0 s ,vr~,0 i ,7r + ,0 s ,7r-,0] , (25) 

where s = (2M - 1) and t = (M 2 -AM- 3). 

The M-indexed family of these cluster states is universal for CV one-way QC. To see this, first 
cut open the toroidal lattice as shown in Figure 8, to form an ordinary lattice. Then, measure q on 
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FIG. 7: Toroidal lattice supergraph and underlying graph structure. Each of the 22 2 macronodes in the 
supergraph (left) consists of four physical nodes, and each color corresponds to one of the four matrix-valued 
weights IP from Eq. (22). Entries in IF specify the real- valued weights connecting the underlying physical 
nodes (right). Measuring q on each physical node in three of the four "layers" leaves the remaining layer in 
a uniformly weighted QC-universal toroidal lattice cluster state. 



three physical nodes per macronode to reduce the supergraph to a uniformly-weighted graph with 
the same structure (see Figure 7), which is known to be a universal graph [18]. These cluster states 
are QC-universal, bicolorable, and orthogonal — but still only block-Hankel. This is nevertheless 
sufficient for simple experimental implementation, using a method we now describe. 

VII. EXPERIMENTAL IMPLEMENTATION 

As already mentioned in Section II, we propose to use optical polarization in order to implement 
the 2x2 blocks tt^. In Eq. (25), each ir^ block corresponds to a single pump frequency. The two 
diagonal elements of in Eq. (19) correspond to interactions between two Y modes and between 
two Z modes. The two off-diagonal elements correspond to YZ and ZY interactions. 

Such couplings can be implemented in a nonlinear crystal such as PPKTP cut along the X axis. 
The point symmetry group of the crystal prescribes which elements of the nonlinear tensor are 
zero [52]. In the case of KTP, the nonzero elements are d2A, c?32) and ^33, which correspond to 
YZY/YYZ, ZYY, and ZZZ interactions respectively. (All other elements involve the X polarization 
and are therefore not accessible in an X-cut crystal.) The difference between tt~ and 7r + is thus 
a 180° phase shift in the Y-polarized pump mode, and a narrowband pump polarized at ±45° in 
the YZ-plane implements a -k^ skew-diagonal band in A. Equation (25) therefore translates into a 
single OPO pumped by exactly 15 frequencies, which must be phase-locked. 

The experimental demonstration of all three interactions simultaneously was performed in a 
PPKTP crystal with a 45.65 /jm period, at a signal wavelength of 1490 nm [33]. The crystal's spatial 
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FIG. 8: Unrolling the torus. The toroidal lattice in Figure 7 can be considered as a particular identification 
of macronodes on an infinite planar lattice that has the local structure of Figure 5. The identification 
is specified by the (nearly) rectangular region delimited by black borders, tessellating the entire plane. 
Macronodes at the same location with respect to their containing region are considered to be the same 
macronode; this enforces the toroidal boundary conditions. Each region has a one-unit "foot" sticking out 
to account for the one- unit twist in each dimension (without the twist, the region would simply be square). 
The region is (M + 1) x (M — 1) + 1 macronodes in size (with the +1 accounting for the "foot"); the 
23 x 21 + 1 region shown corresponds to the 22 2 -macronode lattice in Figure 7. Measuring q on each node of 
the macronodes highlighted red "unrolls" the torus into an ordinary lattice of size M x (M — 2), highlighted 
blue. The M-indexed family of such lattices is universal for QC; thus, the original toroidal lattice is, as well. 



modulation of its nonlinear tensor yielded quasiphasematching of the interactions using different 
harmonics of the same poling period, a remarkable triple coincidence. However, simultaneity at any 
point in the crystal is not required. It suffices to have all interactions occurring every roundtrip 
in the cavity, since the fields undergo many of these before being coupled out, which leads to 
integration of the quantum fluctuations over the OPO cavity's storage time and to the well-known 
squeezing spectrum [53]. (This does yield a finite bandwidth for quantum processing and can 
be circumvented by the use of short-cavity OPOs in order to increase the cold cavity linewidth.) 
Hence, a suitable nonlinear material can be straightforwardly designed by quasiphasematching 
each interaction independently, with its own poling period in separate regions of the crystal. This 
will allow us to use the first (most efficient) modulation harmonic for each interaction, as well as 
to tailor the respective lengths of the poled regions so as to compensate for the different values 
of the aforementioned nonlinear coefficients, and achieve the equal theoretical coupling strengths 
used in this paper. Note that unequal interaction strengths can be compensated by varying the 
relative pump power between the two polarizations, thereby providing a fine tuning mechanism for 
optimizing the entanglement generation. These design improvements to the nonlinear crystal are 
in progress. 

Also, interactions with OFC modes outside the desired subset must also be strictly suppressed 
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in order to realize the exact desired graph state. This can be done by cavity mirror design and/or 
by quasiphasematching bandwidth design [32] and does not present a fundamental issue. While 
the pump spectrum is relatively complex — requiring 15 frequencies — that number is constant with 
respect to the lattice size, making this construction extremely scalable. 

Note that, in principle, additional physical parameters like wave vector direction or transverse 
mode structure could be used to directly implement block-Hankel A's with larger blocks (e.g., 4x4). 
This would reduce pump complexity but require a more sophisticated OPO, in which the nonlinear 
couplings between these additional quantum labels should be precisely defined. Other avenues, such 
as implementation of nonlinear couplings in slow light media, are conceivable. Of particular interest 
would be the use of nonlinear waveguides, seeded by the OFC from a modelocked laser. The interest 
here is to obtain nonlinear efficiency via the intensity enhancement from optical confinement, while 
doing away with the resonant buildup of intracavity losses in an OPO. Recently, an improved, 
though still modest, amount of -2.2 dB squeezing was obtained in a PPKTP waveguide [54]. 

Finally, we examine the experimental requirements for scaling the CV cluster state with N ver- 
tices. As we have shown, the pump spectrum is independent of N and the OFC has constant and 
large size. The scaling requirements are twofold. 

First, the number of pump photons must increase with the number of entangled qumode pairs, 
or graph edges. The toroidal square lattice is a sparse graph, with 0(N) edges for N vertices, 
in lieu of the N(N — l)/2 edges of a complete graph. Increasing the number of vertices therefore 
yields a linear increase of the number of edges, i.e. of the entangled qumode pairs. Thus the 
number of pump photons scales linearly with the number of entangled modes. Since a two-mode 
squeezing experiment in a material like KTP can utilize as low as a few milliwatts of pump power 
and since KTP can withstand up to a few watts of pump power (continuous- wave, focused), this 
allows scaling by 3 orders of magnitude. This figure is only given as an example of what the 
possibilities are, without any optimization, in a typical nonlinear material with a fairly high (but 
not the highest) optical damage threshold. 

Second, the bandwidth of the nonlinear coupling must encompass all entangled qumodes, that 
is, as the graph grows, the skew diagonals in A grow, again linearly with N. The phasematching 
bandwidth is inversely proportional to the length of the nonlinear crystal (over which the nonlinear 
effect is integrated) hence gaining bandwidth by shortening the crystal would yield a reduction of 
the coupling strength and an optimum must be found. Examining again — and conservatively — 
typical figures for existing experiments gives a phasematching bandwidth of 100 GHz to 1 THz to 
be compared to an OPO FSR of 100 MHz (1.5 m long cavity) gives of the order of 10 3 -10 4 modes. It 
is therefore clear from the present estimations that our approach already has a significant realistic 
scaling potential, in line with the theoretical promise. 

VIII. MEASUREMENTS 

Up to this point we have focussed entirely on creating the universal resource for one-way quan- 
tum computation, a CV cluster state. Of course, any useful implementation also needs a viable 
scheme for performing local adaptive measurements. In the original paper on CV one-way quantum 
computing [18], the authors show that universality in an optical measurement-based scheme can 
be achieved with just two types of measurement: (1) homodyne detection and (2) a non-Gaussian 
projective measurement such as photon counting. Homodyne detection is easily performed at close 
to unit efficiency, and the best photon number-resolving detectors currently operate around 95% 
efficiency [55], with the technology continually improving. As a further simplification, recent results 
show that all photon-counting measurements can be performed first, transforming the cluster state 
into a non-Gaussian resource state [56]. This is useful because all remaining measurements are 
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homodyne detections, which are considerably simpler. These measurements are no longer paral- 
lelizable [18, 57] because they are used to gate-teleport non-Gaussian operations onto the encoded 
quantum information in the cluster. 

The unique aspect of our scheme as compared to other optical implementations using CV clus- 
ter states is that our method packs all qumodes into a single beam. Thus, separately addressing 
the polarization- and frequency-labeled qumodes is an obstacle not present in implementations 
relying on spatially separated beams [19, 47, 48, 49]. While separating orthogonal polarizations 
is straightforward, doing so with neighboring cavity modes may seem difficult at first glance since 
an ordinary prism would not have enough resolving power to discriminate optical waves separated 
by radio frequencies. To get around this problem, one can take advantage of new tools for optical 
frequency combs such as the virtually-imaged phased array (VIPA) disperser, which is an interfer- 
ence filter that provides angular separation at the needed resolution and is compatible with lossless 
quantum optical designs [58, 59, 60]. Polarization optics and VIPAs can therefore be used to sep- 
arately address each qumode with individual detectors (whose number necessarily scales linearly 
with the size of the quantum register in any implementation). 

IX. ERROR CORRECTION AND FAULT TOLERANCE 

The main difficulty with any quantum computing scheme is preventing errors and decoherence. 
The theory of quantum error correction and fault tolerance [2] was developed for exactly these 
reasons. While most of the existing literature deals with the case of qubit systems, or more 
generally, finite dimensional systems, our scheme must address the additional difficulties of working 
with continuous variables. Fortunately, error correcting codes have been developed in the CV 
setting [36, 37, 61, 62, 63, 64], usually by mapping the CV system into an effective qubit space [37, 
65]. In some cases, this can be used to establish a fault-tolerance threshold [66]. Lower-dimensional 
encodings may be the only option for fault tolerance using CV cluster states, given the analogy 
with classical analog computation, where no error threshold exists. 

The finite-squeezing approximation is a special case of these more general considerations of 
error correction and fault tolerance. Certainly, more squeezing is preferable to less, but the amount 
required for any particular QC task remains an open question. Consequently, it's unclear how the 
nonlinear coupling strength k in Eq. (2) will need to scale with N for any particular QC application. 
While finite squeezing errors can be mitigated in medium-sized proof-of-principle experiments with 
CV cluster states [18], the task of rigorously establishing a fault-tolerance threshold for CV one-way 
QC is an important open problem. We hope to spur further investigations along these lines. 

X. CONCLUSION 

Quantum computation has come a long way in the last 15 years. There are now multiple 
schemes — both experimental and theoretical — for processing quantum information, all of which 
show promise in different areas. The particular scheme proposed here combines the power of 
cluster-state quantum computation, which harnesses entanglement to substitute for coherent uni- 
tary evolution, with the unorthodox model of quantum computation using continuous variables. 
Experimentally, we propose an efficient method of generating very large continuous-variable cluster 
states using an exceptionally compact experimental setup: a single optical parametric oscillator 
does all the work, entangling all the desired modes at once. One surprising feature of this scheme 
is that the complexity of the required pump spectrum is constant even as the cluster state is scaled 
up to include more modes. Much of the experimental technology required already exists, and 
implementation is underway at the University of Virginia. 



18 



The entangled states produced by this method will also be objects of interest in the study of 
entanglement at mesoscopic scales. Further simplification of the method may be possible using 
additional degrees of freedom, such as spatial modes. This method of entangling — in one fell 
swoop — a large optical frequency comb into a continuous-variable cluster state is the first "top- 
down" approach proposed for one-way quantum computation using optical encodings of information 
in either discrete or continuous variables. While open questions remain about the effects of finite 
squeezing on scalability for particular quantum computing tasks, the unprecedented scalability of 
this method is exciting and motivates further theoretical and experimental research. 
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